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Abstract 

We investigate the behavior, as £ ^ 0, of the nonlocal Allen-Cahn 

equation ut = Am H — -f{u,e / u), where f{u,0) is of the bistable 

type. Given a rather general initial data uq that is independent of £, 
we perform a rigorous analysis of both the generation and the motion 
of interface, and obtain a new estimate for its thickness. More precisely 
we show that the solution develops a steep transition layer within the 
time scale of order ln£|, and that the layer obeys the law of motion 
that coincides with the limit problem within an error margin of order e. 

Key Words: reaction-diffusion equation, nonlocal PDE, singular per- 
turbation, motion by mean curvature 0. 



1 Introduction 

This paper is concerned with the singular limit, as e — > 0, of the nonlocal 
Allen-Cahn equation 



in 



Ut = An H — 2 / ^ y ^) in X (0, oo) 
du 

— = on dQ X (0, oo) 

^u{x , 0) = uo{x) in J7, 



where is a smooth bounded domain in (A^ > 2) and u the Euclidian 
unit normal vector exterior to dQ. We assume that the nonlinearity f{u,v) 
is smooth and that f{u) := f{u, 0) is given by f{u) := —W'{u), where W{u) 



^AMS Subject Classifications: 35K57, 45K05, 35B25, 35R35, 53C44. 
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is a double-well potential with equal well-depth, taking its global minimum 
value at u = ±1. More precisely we assume that / has exactly three zeros 
— 1 < a < 1 such that 

/'(±1) < 0, /'(a) > (bistable nonlinearity) , (1.1) 

and that 

j f{u)du = 0. (1.2) 

The condition (jl.ip implies that the potential W{u) attains its local minima 
at u = ±1, and (|1.2p implies that W{—1) = W{+1). In other words, the 
two stable zeros of / have "balanced" stability. 

Concerning the initial data uq, we assume its smoothness and choose 
Co > 1 such that 

\\^o\\co(n) + l|V'"o|lco(n) + ll^^^o|lc70(n) < Co- (1-3) 
Furthermore we define the "initial interface" Fq by 

Fo := {x e r2| Uo{x) = a}, 

and suppose that Fq is a smooth closed hypersurface without boundary, such 
that, n being the Euclidian unit normal vector exterior to Fq, 

Fo CC n and Vuo{x) / if x G Fo, (1.4) 

Uq > a in Qq, uo < a in , (1-5) 

where Qq denotes the region enclosed by Fq and rj^J" the region enclosed 
between dQ and Fo- 

Before going into more details, let us recall known facts concerning the 
"usual" Allen- Cahn equation, namely 

ut = Au+ ^f{u). 

The singular limit was first studied by Allen and Cahn [1] and by Kawasaki 
and Ohta [H]. By using formal asymptotic arguments, they show that the 
limit problem, as e — > 0, is a free boundary problem: the motion of the limit 
interface is ruled by its mean curvature. More precisely, the solution of 
the Allen-Cahn equation tends to a step function taking the value +1 on one 
side of an moving interface, and —1 on the other side. This sharp interface, 
which we will denote by Fj, obeys the law of motion = — k, where Vn is 
the normal velocity of Ft in the exterior direction and k the mean curvature 
at each point of Ft. 

Then, some rigorous justification of this procedure were obtained. In the 
framework of classical solutions, let us mention the works of Bronsard and 
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Kohn [7], X. Chen [8l[9], and de Mottoni and Schatzman [T8l[T9]. Later, in 
[3], the authors prove an optimal estimate for this convergence for solutions 
with general initial data. By performing an analysis of both the genera- 
tion and the motion of interface, they show that the solution develops a 
steep transition layer within a very short time, and that the layer obeys 
the law of motion that coincides with the formal asymptotic limit Vn = —k 
within an error margin of order e (previously, the best thickness estimate in 
the literature was of order e| lne|, [8]). For similar estimates of the thick- 
ness of the interface in related problems we refer to [T] (reaction-diffusion- 
convection system as a model for chemotaxis with growth), [15] (inhomo- 
geneous Lotka-Volterra competition-diffusion system), [2] (fully anisotropic 
Allen-Cahn equation). 

Since the classical motion by mean curvature may develop singularities 
in finite time (extinction, "pinch off' phenomena...), one has to define a 
generalized motion by mean curvature in order to study the singular limit 
of the Allen-Cahn equation for all time. One represents Tt as the level set of 
an auxiliary function which solves (in the viscosity sense) a nonlinear partial 
differential equation. This direct partial differential equation approach was 
developed by Evans and Spruck [13], Chen, Giga and Goto [H]. In this 
framework of viscosity solutions, we refer to Evans, Soner and Souganidis 
[12], Barles, Soner and Souganidis 0, Barles ans Souganidis [6], Ilmanen 
|16j for the singular limit of reaction-diffusion equations, for all time. 

We now turn back to the nonlocal Allen-Cahn equation. Problem (-P^) 
was considered by Chen, Hilhorst and Logak [TU]. In order to underline its 
relevance in population genetics and nervous transmission, they first show 
that {P^) can be seen as the limit, as a — > and r — > 0, of the FitzHugh- 
Nagumo system 



Then, they study the motion of transition layers for the solutions of (-P^). 
More precisely, for "well-prepared" initial data, they prove that, as e — > 0, 
the sharp interface limit, which we will denote by Ft, obeys the law of motion 



where Vn is the normal velocity of F^ in the exterior direction, k the mean 
curvature at each point of Fj, Q'j' the region enclosed by Fj, the region 
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enclosed between 9fi and r^, cq the constant defined by 



/ —{u,0)du 
Co = —I , (1-6) 

[2{Wiu) -W{-l))]^^^du 



and \A\ the measure of the set A. As explained in [TO], the Problem {P^) 
possesses a unique smooth solution locally in time, say on some [0, T]. More- 
over, in contrast with the "usual" motion by mean curvature which shrinks 
in finite time, the nonlocal effect allows the possibility of nontrivial station- 
ary state (see [T7] for a discussion in the radially symmetric case). 

The goal of the present paper is to make a detailed study of the limiting 
behavior of the solution of Problem (P^), without assuming that the ini- 
tial datum already has a a specific profile with a well-developed transition 
layer. In other words, we study the generation of interface from arbitrary 
initial data. Moreover, we obtain an improved error estimate, of 0(e), be- 
tween the solutions of (P^) and those of (P^)- 

Our main result, Theorem ll.il describes the profile of the solution after 
a very short initial period. It asserts that: given a virtually arbitrary initial 
data uq, the solution u'^ quickly becomes close to ±1, except in a small 
neighborhood of the initial interface Tq, creating a steep transition layer 
around Tq {generation of interface). The time needed to develop such a 
transition layer, which we will denote by t^, is of order lne|. The theorem 
then states that the solution remains close to the step function H on the 
time interval {motion of interface), where u is defined by 

I —1 in ^7 r -, , - 

u{x,t) = { I forte[0,r]. (1.7) 

I -|-1 m il^ 

In other words, the motion of the transition layer is well approximated by 
the limit interface equation {P^)- 

Theorem 1.1 (Generation, motion and thickness of transition layers). Let 

T] be an arbitrary constant satisfying < r/ < min(a -|- 1, 1 — a) and set 

/i = f'{a). 

Then there exist positive constants Eq and C such that, for all e G (0,eo) 
and for all <t <T, where t^ := /^~^e^| lne|, we have 

' [-1-77, +1+77] if x£Mce{rt) 

u'{x,t)£ h-i- 7], -1 + 7]] if xen7\j\fce{rt) (1-8) 
J+1-77,+1 + 77] if xeni\Afce{rt), 

where J\fr{Tt) '■= {x G Q,dist{x,Tt) < r} denotes the r -neighborhood ofVt- 
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The estimate (II. 8p implies that, once a transition layer is formed, its 
thickness remains within order e for the rest of time. 

Corollary 1.2 (Convergence). As e ^ 0, converges to u everywhere in 
Uo<t<T(f^? X {t}). 

This paper is organized as follows. In Section [2] we study the generation 
of transition layers that takes place in a very short time range. Section [3] 
is devoted to the construction of a pair of sub- and super-solutions for the 
study of the motion of interface. In Section H] by fitting the pair of sub- and 
super-solutions of Section [2] into the pair of Section [3l we prove our main 
result. Theorem 11.11 Since our arguments rely on a nonlocal comparison 
principle borrowed from p!0], we recall it in a short appendix. 



2 Generation of interface 

In this section, we investigate the generation of interface, namely the rapid 
formation of internal layers that takes place in a neighborhood of Tq = 
{x G Q,\ uo{x) = a} within the time span of order e^|lne|. In this earlier 
stage, the diffusion term is negligible and the partial differential equation is 

approximated by the nonlocal equation ut = —^f{u,e / u) and so, by the 

Jn 

ordinary differential equation 

ut = \{f{u) + 0{e)). (2.1) 
In the sequel, 7]q will stand for the quantity 

rjo := min(a +1,1 — a). 
The main result of the present section is the following. 

Theorem 2.1 (Generation of interface). Let r] G (0,ryo) be arbitrary and 
define fi as the derivative of f{u) at the unstable zero u = a, that is 

/i = /'(a). (2.2) 

Then there exist positive constants Eq and Mq such that, for all e £ (0,eo), 

(i) for all X E O, 

-l-r?<ti^(x,/i~^e2|lne|) < 1 + r/, (2.3) 

(ii) for all X Q such that \uq{x) — a\ > Mqe, we have that 

if uo{x)>a + MQe then {x, ^i~^e^\\ne\) > I - r], (2.4) 
if uo{x)<a-Moe then u^(x, lne|) < -1 + (2.5) 
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The above theorem wih be proved by constructing a suitable pair of 
sub- and super-solutions based upon the ordinary differential equation (j2.ip . 
Note that the assumption of balanced nonlinearity ()1.2p is useless for the 
proof of the generation of interface property. 



2.1 The bistable ordinary differential equation 

We first consider a slightly perturbed nonlinearity: 

Mn) ■.= f{u) + 6, 

where 5 is any constant. For \5\ small enough, this function is still of the 
bistable type. More precisely, if Jq is small enough, then for any 6 S 
(— (5o,5o), fs has exactly three zeros, namely a-{5) < a{6) < a+{5), and 
there exists a positive constant C such that 

|a_((5) + 1| + \a{6) - a\ + \a+{d) - 1| < C\6\, (2.6) 

H6)-fi\<C\6\, (2.7) 

where 

fi{6) := f's(a{5)) = f'{a{5)). 

Now for each 5 G {—5o,6o), we define Y{t,S,;5) as the solution of the 
ordinary differential equation 

iYr{T,C;5) =fs{Y{T,C;5)) for r > 

\ y(o,e;5) =C, 

where ^ varies in (— 2Co,2Co), with Cq being the constant defined in (jl.3p . 
We claim that Y{t,£,;6) has the following properties. 

Lemma 2.2. There exist positive constants 5q and C such that, for all 
(T,e,<5) G (0,oo) X [-2Co,2Co] x [-60,60], 

(i) \Y{r,C;6)\<2Co 

(ii) 0<Y^{t,C-6) 

(iii) |^(r,e;<5)|<C(e'^(^)--l). 

Property (i) is a direct consequence of the profile of fs and so of the 
qualitative properties of the solution of the bistable ordinary differential 
equation (|2.8p : for proofs of (ii) and (iii) we refer to [3], subsection 4.1. □ 
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2.2 Construction of sub- and super-solutions 

We are now ready to construct a pair of sub- and super-solutions in order 
to prove the generation of interface property. By using some cut-off initial 
data (see [3], subsection 3.2) we can modify slightly uq near the boundary 
90 and make, without loss of generality, the additional assumption 

^ = on dn. (2.9) 

We set ^ ^ 

wf{x,t) = y(^^,uo{x) ±e^r{±eg, ^);±e^^, 

where the function r(5, r) is given by 

r((5,r) =a(e^('^)^-l), 

and the constant G by 

df 

Q = 2Co\^\ max \ — (u,v)\. 

{u,v)G[-2Co,2Co]x[-l,l] OV 

Lemma 2.3. There exist positive constants Eq and such that, for all 
£ G (0, So); {"^7 ^''^t) ^'^ P^'^''^ of sub- and super- solutions for Problem (P^), 
in the domain Q x (0, lne|). 

Before proving the lemma, we remark that w~{x, 0) = wf{x, 0) = uq{x). 
Consequently, by the comparison principle, we obtain 

w-{x,t) <u^x,t) <wf{x,t) for ah n X [0,/i-^e2|lne|]. (2.10) 

Proof. First, the inequality w~ < wf follows from the fact that Y{t, ^; 6) 
increases with both (see (ii) Lemma 12. 2p and 6 (as easily seen from the 
ordinary differential equation). Next, (j2.9p implies that both and w~ 
satisfy the Neumann homogeneous boundary conditions. Hence, it remains 
to prove the inequalities C^wf > and C^w~ < (see Definition lA.ip . 
provided that the constants eo and are appropriately chosen. 

If eo is sufficiently small, we note that zizeQ £ {—5o,5o) and that, in the 
range < t < //"""^e^l lne|, 

|e2c^(eA'(±^5)*/=' - 1)1 < e2(7^(£-M(±eG)/M _ i) < Cq, 
using (12. 7p . The above inequality implies 

uoix)±e\{±£g,-^) £ [-2Co,2Co]. 
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These observations allow us to use the results of Lemma [2.21 with the choices 
r := t/e^, ^ := uo{x) ± e^r(ibeC/, t/e^) and 6 := ±eQ. In particular, it fol- 
lows from property (i) that | / w^{x,t)dx\ < 2Co\Q\ which in turn implies 

Jn 

(thanks to the choice of G) 



_max f{wf,£s)<f{wf) + eg. 



(2.11) 



In view of the above inequality, some straightforward calculations yield 

where the argument ^—,^0(2^) + e^Cj,(e'^*-^^^*''^^ — X)\eQ^ of the function 

Y and its derivatives is omitted. Noticing that the ordinary differential 
equation (|2.8p writes as Y^ = f(Y) + eQ, we get 



C+wt > n 



Ci,ii{eg)e 



Auo 



Y, 



\Vuo\' 



Using the estimate (iii) in Lemma 12.21 we obtain 

£+w+ > Y^ fa//(eg)e'^(=^)*/^' - \Auo\ - C(e'^(^^)*/=' - l)|Vuo|2 



> Y, 



iCMeg) - C\Vuo\')e 



|Ano| +C|Vno|^ 



In view of ()2.7p . this inequality implies that, for e G (0,eo)) with small 
enough. 



C+wt > Yf 



C*-fJ. — CCq^ — Co 



>o, 



by choosing large enough. 

Since one can prove C^w~ < by similar arguments, this completes the 
proof of Lemma 12.31 □ 



2.3 Proof of the generation of interface property 

In order to prove Theorem 12.11 we first quote a lemma from [3] ; it makes 
more precise the bistable behavior of the ordinary differential equation by 
giving basic estimates of the function Y{t,£^; ^eQ) at time r = lne|. 

Lemma 2.4. Let r] € (0,?7o) be arbitrary; there exist positive constants Eq 
and Mi, such that, for all e G (0, eo)? 



(i) for aine{-2Co,2Co), 

- 1 - ri <Y{fi~^\ In e\,^;±£g) < 1 + ??; 



f2.12) 
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(ii) for all ^ G (— 2Co, 2Co) such that |^ — a| > M^^e, we have that 

if i>a + M^e then lne|, ^; ieC?) > 1 - r/, (2.13) 

if (<a-M^e then lne|, ^; ±eC?) < -1 + ??• (2.14) 



Proof of Theorem 12.11 By setting t = fi lne| in ()2.10p . we get 

Y(^fi^^\ lne\,uo{x) - e^r{-eg, fj.~^\ lne|); -eG 

< u%x,fi~'^e^\lne\) < Y (n~^\lne\,uo{x) + e'^r{eG, n~'^\lne\);+eg 



(2.15) 

We note that (12. 7p imphes 

^.^M^M±££)ine = 0, (2.16) 

so that, if £o is sufficiently smaU, 

e\{±eg,fi-'\lne\) = ae(e('^-^(±^^»/^ - e) e (^C.e, ^C,e), 

and, for all x € Q, it holds that uo{x) ± £'^r{±£Q , fM~'^\ In e\) £ (-2Co, 2Co). 
Hence, the result (j2.3p of Theorem 12. II is a direct consequence of (|2.12p and 
(I2J5]) . 

Next we prove (|2.4p . We take x £ il. such that uq{x) > a + Mqe; then 

tio(x) - e2r(-eg,^-i| lne|) > a + Moe - fC^e 

> a + M^e, 

if we choose Mq large enough. Using ()2.15p and ()2.13p we see that inequal- 
ity (j2.4p is true. The inequality (j2.5p can be shown the same way. This 
completes the proof of Theorem 12.11 □ 

3 Motion of interface 

In Section^ we have proved that the solution of Problem (P^) develops a 
clear transition layer within a very short time. The aim of the present section 
is to show that, once such a clear transition layer is formed, it persists for the 
rest of time and that its law of motion is well approximated by the interface 
equation (P^)- In order to study this latter time range where the motion of 
interface occurs, we will construct another pair of sub- and super-solutions 
{u~,uf) for Problem (-P^). To begin with we present mathematical tools 
which are essential for this construction. 
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3.1 Preliminaries 



The "cut-ofF signed distance function". Let T = Uo<t<r(rt x {t}) be 
the solution of the limit geometric motion problem {P^) and let d be the 
signed distance function to F defined by: 

~ (— dist(x,rt) for X G rir 

d{x,t) = \ * (3.1) 

[ dist(x,rt) for X G ri^, 

where dist(3;,rj) is the distance from x to the hypersurface Tt in The 
"cut-off signed distance function" d is defined as fohows. First, choose do > 
small enough so that the signed distance function d defined in (j3.1|) is smooth 
in the following tubular neighborhood of F: 

{{x,t) G X [0,r] I \d{x,t)\ < 3do}, 

and that 

dist(Fi, on) > 3do for ah t G [0, T]. (3.2) 
Next let C(s) be a smooth increasing function on M such that 

s if \s\ < do 

C{s) = < -2do if s < -2do 
2do if s > 2do. 

We then define the cut-ofi' signed distance function d by 

d{x,t) =c(^d{x,t)y (3.3) 



Note that, in view of (j3.2p and the definition of d, the equality Vd = 
holds in a neighborhood of dO,. Note also that the equality |V(i| = 1 holds 
in the region {{x,t) G 17 x [0,T] | |d(a;,t)| < do}. Moreover, since Vd 
coincides with the outward normal unit vector to the hypersurface Ff, we 
have dt{x,t) = —Vn, where Vn is the normal velocity of the interface Tt 
in the exterior direction. It is also known that the mean curvature k of 
the interface is equal to Ad. Hence, since the moving interface F satisfies 
Problem (P"), an alternative equation for F is given by 

dt = Ad- co7(t) on Tt, (3.4) 

where 7(t) := — 

The one dimensional standing wave Uo- Let Uo{z) be the unique solu- 
tion of the stationary problem 

f uo" + fm = 

\Uoi-oo) = -l, [/o(0) = o, C/o(+oo) = +l. 
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This solution represents the first approximation of the profile of a transition 
layer around the interface observed in the stretched coordinates; it naturally 
arises when performing a formal asymptotic expansion of the solution 
(see [3], Section 2). Note that the "balanced stability assumption", i.e. 

the integral condition / du = 0, guarantees the existence of such a 

standing wave. In the simple case where f{u) = u{l — u^), we know that 
Uo{z) = tanh(z/\/2). In the general case, the following standard estimates 
hold. 

Lemma 3.1. There exist positive constants C and A such that 

< 1 - Uo{z) < Ce-^\^\ forz>0 
< Uo{z) + 1 < C7e-^l^l for z<0. 
In addition, Uq is a strictly increasing function and, for j = 1,2, 

\DWo{z)\ < Ce-^l^l for zeR. (3.6) 

The solution Ui of a linearized problem. Let Ui{z,t) be the solution 
of the problem 



where 



Uu, + f'iUoiz))Ui = 7(i) (-|^(t/o(z),0) - coUo'iz) 
C/i(0,t)=0, Ui{-,t)£L 

Of 



(3.7) 



{Uo{z),0)Uo'iz)dz 

c,:=-^-^. . (3.8) 

/ {Uo'f{z)dz 

Again, the above problem arises when performing a formal asymptotic ex- 
pansion of the solution u^. Since (|3.7|) can be seen as a linearized problem for 
p.Sp . its solvability follows from a Fredholm alternative: thanks to the def- 
inition of Co, Uq' turns out to be orthogonal to the right-hand side member 
of (13. 7p . Moreover, there exist constants M > and C > such that 

\Ui{z,t)\<M, (3.9) 

\Uu{z,t)\<C, (3.10) 

\UUz,t)\ + \Uu,{z,t)\ < Ce-^l^l, (3.11) 

for all (z,t) £ R X [0, T]. We omit the details and refer the reader to [3], 
Section 2. 

Note that, by multiplying equation (13. 5p by Uq' and integrating from 
-oo to z, we obtain Uo'{z) = [2{W{Uo{z)) - Wi-l))]^/"^. Using this, it is 
now a matter of routine to deduce from ()3.8p the more intrinsic expression 

dLBD. 
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3.2 Construction of sub- and super-solutions 

We look for a pair of sub- and super-solutions uf for {P'^) of the form 

+ , , f d(x,t) ±ep(t)\ f d(x,t) ±ep(t) \ , 
uf{x,t) = Uo ^ )+eUii ^,t)±q{t), (3.12) 



where 



p(t) = -e-^*/^ + e^* + K, 
q{t) =a{Pe-l^'/'' +e^Le^'). 



Note that q = ae^ pt- It is clear from the definition of that 

lim„f(..t) = (+' f"'J'(-.*)^U«S.S-W>'W) (3.13, 
\ -1 forall(i,t)eUo<t<T(!!rx{*})- 

The main result of this section is the following: 

Lemma 3.2. Choose /3 > and a > appropriately. Then for any K > 1, 
there exist positive constants Oind L such that, for any e G (0, Eq); (^7 ' '^t) 
is a pair of sub- and super- solutions for (P^) in the domain Q x [0,T]. 

Proof. First, we claim that (|A.14p and (|A.16p hold as a consequence of 
()4.5p . Then, since d is constant in a neighborhood of dQ, both uf and u~ 
satisfy the Neumann homogeneous boundary conditions. Hence it remains 
to prove the inequalities C+uf > and C-U'^ < 0, provided that the various 
constants are appropriately chosen. 

We start with some useful inequalities. On the one hand, by assumption 
(jl.ip . there exist positive constants b, m such that 

f'{Uo{z)) < -m if Uo{z) G [-1, -1 + 6] u [1 - b, 1]. (3.14) 

On the other hand, since the region {z G M| Uo{z) G [—1 + 6, 1 — 6] } is 
compact and since Uq' > on M, there exists a constant ai > such that 

Uo'{z)>ai if C/o(z) G [-1 + 6, 1-6]. (3.15) 

We set 

P=j, (3.16) 

and choose a that satisfies 

< cj < min ((Tq, 0-2), (3.17) 

where 

ai 1 4/3 

m + Fi p + 1 f/ (/3 + 1) 
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the constants -Fi and H being given by 

^1 := II/'IIl°°(-i,i), H := ^ max \\Hess(^u,v)f\\, (3-18) 

where ||^|| : = \. Combining (j3.14p and (jS.lSp . and considering 

that (T < (To, we obtain 

Uo'{z) - af'{Uo{z)) >am for - oo < z < oo. (3.19) 



Now let X > 1 be arbitrary. In what follows we will show that C-^-uf > 
provided that the constants and L are appropriately chosen. We recall 
that — 1 < C/q < 1 and that \Ui\ < M. We go on under the following 
assumption 

EqM < 1, elLe^'^ < 1 . (3.20) 
Then, given any e G (0,eo)! since a < ai we have < q{t) < 1, so that 

-3<uf{x,t) <+3. (3.21) 

In order to evaluate the "nonlocal part" of C+uf , we need bounds for 
the quantities / u^{x,t) dx. For the sake of clarity, the arguments of most 

of the functions are omitted in the following. We write / uf dx = / {Uq — 

Jn Jn+ 

1) dx + \nf\+ / {Uo + l)dx - + / (eUi + q) dx, which yields 
Jn^ Jn 

/ ufdx-j{t) =\n\q{t)+ / eUidx 
Jn Jn 

+ [ {Uo-l)dx+ [ {Uo + 1) dx 
Jnf Jn^ 

=: \n\qit) + Ie{t)+I+{t) + I_it). 

In the following we will denote by C various positive constants that are 
independent of e E (0, eo). Since Ui is bounded, we have \Ieit)\ < Ce. 
In order to estimate I+{t), we use the partition 

nf = {x\ d{x,t) > do} U {x| < d{x,t) < do}. 

First assume d{x,t) > do- From Lemma l3. 11 we deduce that 
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from which we infer that 

< / f 1 - f;„(^(M±M!) dx < Ce-^^o/e. 

Jd(x,t)>do \ £ / 

In order to estimate the integral on {x\ < d{x,t) < do}, we use arguments 
similar to those used in [TU]. We denote by J(s,d) the Jacobi of the trans- 
formation X ^ {s,d), where s{x,t) is the projection of x on Tt along the 
normal of Tt and d{x,t){= d{x,t)) is the signed distance defined above; we 
define Cj := maxo<t<T \\J{-,-)\\L°-{rtx[-do,do])- ^^^^ ^^^^"^^ 

< /■ f 1 - f,„(*2^)±£r*>,) <i, 

Jo<d(x,t)<do \ ^ J 



Ft JO 



<Cje/ / (l-C/o(z + p(t))) 



Tt JO 



r- r+oo 

<Cje / (1 - Uo{u)) duds < Ce. 
Jvt Jo 



Tt Jo 

As far as I-{t) is concerned, we first assume that d{x,t) < —do. Note 
that 0<i^ — l<p< e^'^ + K. Consequently, if we assume 

e'^^ + K<^, (3.22) 
2eo 

then — — IpI > By using similar arguments as the ones above, we obtain 
< / ( f7„(^(^2^)±M^) + i)dx< Ce-^^"/(^^). 



'{x,t)<—do 

Concerning the region {x\ — do < d{x,t) < 0}, we get 
< / (u/i^,t)+epit)^ ^ \ 

J -da<d{x,t)<0 \ ^ J 

= £ (f/o(^^^) + l) J{s,r)drds 

< Cje [ r {Uo{z + p{t)) + 1)) dzds 

JVt J~do/e 



r /•do/{2eo) 

Cje / / {Uo{u) + 1) duds < Ce. 

JTt J -oo 



< 
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Since one would obtain similar estimates with replaced by , the 
above estimates yield 



dx — 7(t) 



< Ce + Cq{t), 



which, in turn, implies 



_max f{ut,es) < f{uf,e-f{t)) + Ce^ + Ceq(t) 



and 



(3.23) 



(3.24) 



_min f{u~,es) > f{u~,e-f{t)) - Ce^ - Ceq{t). (3.25) 



Now, we can turn back to the proof of Cj^u'^ > 0. From the above 
inequality, we get 

£+n+ > {ut)t - An+ - l/(n+,e7(0) -C- C-q{t). (3.26) 



Straightforward computations yield 

{u+)t = Uo'ij +pt) + eUit + UiM + ept) + qt 

where the function Uq, as well as its derivatives, are evaluated at 2; = 
[d(x,t) + ep{t)'j/e, whereas the function Ui, as well as its derivatives, are 

evaluated at (^{d{x, t) + ep{t)^ /e, • We also have 

df 

fiut,e^it)) <f{Uo) + (eC/i +g)/'(C/o) +e7(i)^(f^o,0) 

+ H Q(eC/i + qf + ^eV(i) + {eUi + q)e-f{t) 



where H was defined in (|3.18|) . Combining the above expressions with the 
equations (13. 5p and (13. 7p for C/q and Ui, we obtain 



where: 



E2=(% + ^]{l-\Vd\'^ 



\q (f'iUo) + lHq\+ Uo'pt + qt 
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^3 = f — + ] (dt -Ad + co7(t)) 



e 

= eUupt + \ {-Hj{t) - HUi - C) 

E5 = -co7{t)Uiz - \HUi'' - lHj\t) - H^{t)Ui - C 
Eq = eUit- 

In the sequel, we estimate the terms Ei — Eq and denote by Cj various 
positive constants that are independent of e. 

3.2.1 The term Ei 

Direct computation gives 

El = ^ e-^'/'\l - a(3) + Le^\l + e^aL), 

where 

In virtue of (j3.19p . we have 



2 



2 r ^LT\ 



I>am- —H{l3 + e^Le 

Combining this, (j3.16p . (j3.20p and the inequahty a < a2, we obtain / > 2cj/3. 
Consequently, we have 

3.2.2 The term E2 

First, in the region where \d\ < do, we have \Vd\ = 1, hence E2 = 0. Next 
we consider the region where \d\ > do- We deduce from Lemma [3. II and from 
(ISTTT) that 

In view of (j3.22p we have < p < — so that 

2^ 



\E2\ < '^SL^~^<inl{2e) < 
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3.2.3 The term E^, 

Recall that 

(dt - Ad)(2;,t) + co7(t) = on Vt = {x £ Q., d{x,t) = {)}. 

Since the interface Tt is smooth, both Ad and dt are Lipschitz continuous 
near Tt- It follows from the mean value theorem applied on both sides of Tt 
that there exists a constant > such that: 

\{dt - Ad)(x,t) + co7(i)| < N\d{x,t)\ for all (x,t) G J7 x (0,r). 

Applying Lemma 13.11 and the estimate (j3.1ip we deduce that 

|^3| < 2iVC^e-^l^/^+Pl 
e 

< 2iVCmax5eK|e|e-^l«+Pl 

< 2iVCmax(|p|,^). 

A 

Thus, recalling that \p\ < e^^ + if, we obtain 

|i?3| <C3(e^*+i^)+C3', 

where C3 := 2NC and C3' := 2NC/X. 

3.2.4 The terms S4, E5 and Eq 

Since Ui, Uiz, Uu and 7 are bounded, it is a matter of routine to see that 
\Ei\ < CA-Pe-^'/'' +£Le^'), {E^l < C5, {EqI < eCq. 

3.2.5 Completion of the proof 

Collecting all these estimates gives 

£+u+ > - ^)e-/3*/^' + (2a/?L - C73 - eC^L)e^' - Cj, (3.27) 
where Cj := C2 + KC3 + C3' + C5 + Now we set 

which, for eq small enough, validates assumptions ()3.20p and ()3.22p . For eq 
small enough, the first term of the right-hand side of (|3.27p is positive, and 

£+u+ > [apL - C3]e^* - C7 > ^(^PL - C7 > 0. 

The proof of ()A.12p is now complete, with the choice of the constants /?, a 
as in ()3.16p . (I3.17p . Since one can prove (1A.13P by similar arguments, this 
completes the proof of Lemma 13. 2i □ 
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4 Proof of Theorem 11.11 

Let 7] £ (0,r/o) be arbitrary. Choose (3 and a that satisfy (j3.16p . (jS.lTp and 

ap < |. (4.1) 

By the generation of interface property, Theorem 12.11 there exist positive 
constants Eq and Mq such that (12. 3|) , ()2.4p and (12. 5p hold with the constant 
r] replaced by cr/3/2. Since Vuq / everywhere on Tq = {x £ r2| uq{x) = a} 
and since Tq is a compact hypersurface, we can find a positive constant Mi 
such that 

if do{x) > Mi£ then uq{x) > a + Mqe ^^^^ 
if dQ{x) < —Mi£ then ^0(2:) < a — Mqe. 

Here dQ{x) := d{x,Q) denotes the signed distance function associated with 
the hypersurface Tq. Now we define functions H~^{x),H~(x) by 



H+{x) 
H-{x) 



+1 + 0-/3/2 if do{x) > -Mie 

-l + ap/2 if do{x) < -Mis, 

+1 - a(3/2 if doix) > Mie 

-1 - ap/2 if do{x) < Mie. 



Then from ([23]), i^^, ([23]) (with 7] replaced by a/3/2) and (|121), we see 
that 

H^{x) < u%x,fi^^e^\ In e\) < H+{x) for x G 0. (4.3) 
Next we fix a sufficiently large constant K > 1 such that 

[/o(-Mi + K) > 1 - ^ and Uo{Mi-K)<-l + ^. (4.4) 

For this K, we choose Eq and L as in Lemma |3.2[ We claim that 

u-{x,0) < H-{x), i7+(x) < ti+(x,0) for X GO, (4.5) 

with {u~,uf) the pair of sub- and super-solutions defined in (|3.12p for the 
study of the motion of interface. We shall only prove the former inequality, 
as the proof of the latter is virtually the same. Then it amounts to showing 
that 

u~{x,0) = Uo{^^-K)+eUi{^^-K,0)-a{(3+e^L) < H-{x). (4.6) 

Recall that \Ui\ < M. Therefore, by choosing eq small enough so that 
EqM < a (3/ 6, we see that 

ujix,0)< Uo{^-K)-la(3. 

E 
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In the range where dQ{x) < Mie, the fact that Uq is an increasing function 
and the second inequahty in (j4.4p imply 

£ b 2 

On the other hand, in the range where do{x) > Mie, we have 

£ D 

This proves (|4.6|) . hence (|4.5|) is estabhshed. 
Combining ()4.3p and ()4.5p . we obtain 

uj{x,0) < u^{x,fi~^£^\ln£\) < uf{x,0). 

Since (u~,uf) is a pair of sub- and super-solutions for Problem the 
comparison principle yields 

u~{x,t) <u%x,t + f) <uf{x,t) forO<^<^-^^ (4.7) 

where = /i~^e^| lne|. Note that, in view of ()3.13p . this is enough to prove 
Corollary II. 2i Now let C be a positive constant such that 

f/o(C-e^^-ir) > 1- I and Uo{-C + e^^ + K) < -1 + ^. (4.8) 

One then easily checks, using (|4.7p and (j4.ip . that, for £q small enough, for 
< t < T — t*", we have 

if d(x,t)>C£ then u%x,t + 1^) > I - n , ^ 

^ ' ^ - ^ ' > - ' (4.9 

if d{x,t)<-e£ then u=(x, t + t^) < -1 + r?, 

and 

u^(2;,t + t^) G [-1 - + 
which completes the proof of Theorem II. 1[ □ 



Appendix - Comparison principle 

The definition of sub- and super-solutions is the one proposed by Chen, 
Hilhorst and Logak [10]. It involves simultaneously a super- and a sub- 
solution. 

Definition A.l. Let {u~ ,u^) be a pair of smooth functions defined on 
Cl X [0, T] and satisfying 

< u+ in X [0,r], (A. 10) 
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and 

^<0<^ on 9J^x(0,T). (A.U) 

We say that n^) is a pair of sub- and super-solutions for Problem {P^) 
if 

:= (n+)t - Au+ - ^ _max /(n+, es) > in 17 x (0, T), 

^ (A.12) 
C-U~ := {u~)t — Au~ 2 — in n X (0,T). 

(A.13) 

As proved in [TU], the following comparison principle holds. 

Proposition A. 2. Let a pair of sub- and super-solutions be given. Assume 
that, for all x £ il, 

u~{x,0) < no(x) < ti+(x,0). (A.14) 

Then, if we denote by the solution of Problem {P^), the function 
satisfies 

uj{x,t) <u^{x,t) <uf{x,t), (A. 15) 

for all {x,t) G X [0,r]. 



As easily seen from the proof in [10], one could replace the assumption 
(|A.10p by the assumption (jA.14p together with the condition that 

/ u^{x,t)dx< / uf{x,t)dx, (A.16) 
Jn Jn 

for all t G [0, to] with to > 0. More precisely if fOill . (fXlG]) . fOT]l . i^KJ2\i 
and (jA.lSp hold, then the conclusion (|A.15p follows. 
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